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ABSTRACT 


Spacecraft operating in low earth orbit face a significant orbital debris impact hazard. 
Of particular concern, in the case of the Space Shuttle, are impacts on critical components 
of the thermal protection system. Recent research has formulated a new material model 
of reinforced carbon-carbon, for use in the analysis of hypervelocity impact effects on the 
Space Shuttle wing leading edge. The material model has been validated in simulations 
of published impact experiments and applied to model orbital debris impacts at velocities 
beyond the range of current experimental methods. The results suggest that momentum 
scaling may be used to extrapolate the available experimental data base, in order to predict 
the size of wing leading edge perforations at impact velocities as high as 13 krn/s. 

NOMENCLATURE 


d, shear damage variable 
D , projectile diameter 
D p , perforation diameter 
D c: coating spall diameter 

D. rate of deformation tensor 

e, Euler parameter vector 

E. deviatoric strain tensor 
E e , elastic strain tensor 
E p , plastic strain tensor 
e p , effective plastic strain 
e, deviatoric strain rate 

F. deformation gradient tensor 


Sij, Kronecker delta 
L. velocity gradient tensor 
N e , number of elements 

R, rotation matrix 

S, deviatoric stress tensor 
S p , effective stress tensor 
s, entropy density 

v. impact velocity 
Y, yield stress 
•0, strain energy density 
9, temperature 
< i i>, impact obliquity 
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INTRODUCTION 


Carbon-carbon composites offer an unusual combination of thermal and mechanical prop- 
erties. 1 Their light weight and high temperature strength satisfy some very stringent design 
requirements for reusable orbital vehicles. 2 The wing leading edge of the Space Shuttle, 
subject to severe thermal re-entry loads, is constructed of reinforced carbon-carbon (RCC) 
panels, coated in silicon carbide to prevent oxidation. 3 Although the thermal properties of 
RCC composites are well understood, 4 much less is known about their dynamic mechanical 
properties. The loss of the Space Shuttle Columbia, 5 apparently due to impact damage on 
the wing leading edge, has motivated recent experimental 6 and computational' work aimed 
at developing a better understanding of the impact response of thermal protection materials. 

The wing leading edge damage to Columbia was unexpected, the result of a relatively 
low velocity impact by a relatively low density projectile. 8 Another impact damage hazard, 
due to space debris in low earth orbit, has long been recognized. This threat involves 
projectiles of very low mass, but much higher density, and impact velocities as high 15 km/s. 
The debris shielding on the International Space Station is designed to defeat centimeter 
sized aluminum projectiles. Although the likelihood of such a projectile striking the Space 
Shuttle is quite low, orbital debris damage by much smaller projectiles is routinely observed 
during post-mission inspections of the vehicle. As a result previous experimental research 
has investigated the response of Space Shuttle thermal protection materials to orbital debris 
impact by spherical aluminum projectiles as large as 0.628 cm in diameter. 9 

Due to the high cost of carbon-carbon composites and the long fabrication lead times 
associated with the preparation of test samples, impact testing of RCC materials has been 
limited. In addition, the limitations of current experimental technology preclude hyperveloc- 
ity impact testing over the entire projectile mass and kinetic energy range of interest. As a 
result, numerical simulation can serve as an important complement to experimental studies 
of the impact response of RCC materials. Numerical models validated by comparison with 
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experiment at velocities below 8 krn/s can be used to extrapolate results into a higher veloc- 
ity impact regime. A coordinated experimental and computational approach to the study of 
RCC response to insulating foam impacts has proven to be productive; 10 the present paper 
extends the latter computational work, to projectiles and impact velocities associated with 
orbital debris impact. In particular it develops a new anisotropic, rate-dependent material 
model for reinforced carbon-carbon, validates that model in three dimensional simulations of 
published hypervelocity impact experiments, and applies the validated formulation in simu- 
lations of impacts at velocities beyond the experimental range. The results indicate that a 
momentum scaling approach used to correlate the available experimental impact data may 
be extrapolated to describe RCC perforation by hypervelocity projectiles at velocities as 
high as 13 krn/s. 

The present paper is organized as follows. The first section outlines the hybrid particle- 
finite element method used in the present study, including the imbedded large deformation 
kinematics and general functional forms for the associated constitutive relations. The second 
section discusses published experimental results on the properties of RCC. The third section 
develops an RCC constitutive model, formulated for use in hypervelocity impact applica- 
tions and reflecting important mechanical characteristics described in the material testing 
literature. The fourth section validates and applies the developed model in a series of three 
dimensional impact simulations. The last section presents conclusions and suggestions for 
related future work. 


NUMERICAL METHOD 


The material model described in this paper was developed for application in a specific 
numerical framework, the hybrid particle-finite element formulation of references 11 and 
12. In order to provide appropriate context, this section summarizes the latter numerical 
formulation, details certain element level kinematics, and provides functional forms for the 
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required constitutive relations. The kinematic and constitutive modeling framework assumed 
here has wide scope, so that the material model described in the present paper may be 
adapted for use in shock physics codes which are based on alternative numerical modeling 
schemes . 13 

The hybrid particle-finite element model employed here takes an explicit state space form. 
The state equations consist of evolution equations for the following variables: 

• translational and rotational momentum vectors for the three dimensional motion of 
ellipsoidal particles, 

• center of mass position vectors and Euler parameters for the particles, the latter pro- 
viding a singularity free description of particle rotations, 

• density and entropy for each particle, and 

• damage and plastic internal state variables for each finite element. 

The state equations are derived using a thermomechanical formulation of the Lagrange equa- 
tions. All inertia effects are modeled using the particles, whose mass centers are also nodal 
coordinates for the finite elements. The volumetric thermomechanical response of the mod- 
eled medium is described by an equation of state for the particles, which may take either an 
analytic or tabular form. 

The material modeling work described in the present paper develops two specific compo- 
nents of the general numerical formulation: 

• a strain energy density in shear, one part of the thermomechanical Lagrangian for the 
modeled particle-element system, and 

• a plasticity model which specifies evolution equations for the plastic internal state 
variables, equations which serve as nonholonomic constraints on the system level model. 
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The strain energy density in shear takes the general functional form 


F = if>(d,e,E,W) (1) 

where d is a shear damage variable, E is the total deviatoric strain, E p is the plastic strain, 
and e is a vector of Euler parameters which relates a material reference frame for each 
element to a single global Cartesian reference frame. The evolution equations for the plastic 
strain components take the general functional form 

E p = E p (s, d, e p , e, J, e, E. E p ) (2) 

where s is an entropy density, e p is the effective plastic strain, e is a deviatoric strain rate, 
and 

J = det{ F) (3) 

where F is the deformation gradient tensor. 

The strain and strain rate variables which appear in the preceding functional forms are 
defined by the following large deformation kinematics. 14 The deviatoric strain is 

E=t(C-l) (4) 

where 

C = F i F. F = (detF)-3F (5) 

The elastic shear strain is defined as 

E e = E - E p (6) 

where the flow rule for the plastic stain tensor must satisfy the isochoric plastic deformation 
constraint 

tr ( C P ~ T & ) =0, C p = I + 2 E p (7) 
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The effective plastic strain is determined by integrating the rate relation 


i p = 



( 8 ) 


with the indicated invariant operator defined by 

1 


IT 1 1 = 


- tr (T T T) 


1/2 


(9) 


for any second order tensor T. The deviatoric strain rate is 


e = | |D'| |, 

where D is the rate of deformation tensor 

1 


D ; = D — -tr( D) I 


( 10 ) 


D 




L = FF 


-l 


( 11 ) 


with L the velocity gradient tensor. 

In the case of anisotropic materials, the constitutive response is described in a material 
reference frame. Here an Euler parameter vector 


e = [e 0 ei e 2 e 3 ] r , e T e = 1 (12) 

is used to define a rotation matrix (R) for each element 



R 
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which relates a material reference frame in each element to the global Cartesian system used 
in the numerical simulations. The rotation matrix relates vector components p in the global 
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coordinate system to corresponding components q described in the material reference frame, 
using 

p = R q (16) 

The corresponding transformation relation for second order tensors is 

P = RQR t (17) 

where the components P refer to the global frame and the components Q refer to the material 
reference frame. 


REINFORCED CARBON-CARBON 


The published material property data base for carbon-carbon composites is limited by 
material costs and proprietary considerations. On the other hand, the complex nature of 
both the material and the application of interest here means that a rather wide range of 
experiments are needed to fully characterize its constitutive response. This section discusses 
some properties of RCC of particular significance in hypervelocity impact applications. 

The most directly relevant experimental results are those of Lu et ah, 6 who performed 
tests at Sandia National Laboratories on samples taken from Space Shuttle wing leading edge 
panels, in support of the Columbia accident investigation. They provide data on elastic mod- 
uli as well as strength measurements obtained from tension, bending, and compression tests. 
Although the elastic moduli measured in tension and compression were similar, strength in 
compression was approximately double that in tension. In addition they reported a strain 
rate dependence of the tensile strength, observing a fifteen percent increase in strength as 
the loading rate increased from 1 to 200 sec -1 . Finally they noted that removal of the silicon 
carbide coating from the tested samples showed little effect on the measured mechanical 
properties. 
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Several different authors have reported results of shear tests performed on carbon-carbon 
composites . 15-20 In the case of the RCC, interlaminar shear strength and stiffness is of 
interest, since oblique hypervelocity impacts will in general lead to multi-axial loading. The 
published shear test data show that interlaminar shear stiffness and strength can differ by 
factors of approximately two and four respectively, from their in-plane counterparts. 

Perhaps the most unusual property of RCC is its increase in strength with tempera- 
ture , 17,21 by as much as a factor of two, as compared to the thermal softening response 
observed in metals. The high temperature strength of RCC is important in hypervelocity 
impact applications, due to the adiabatic heating typically associated with shock loading. 

The preceding references, along with the equation of state literature 22 and published data 
on the thermal properties of RCC , 4,23 were used to estimate the material parameters used 
in the simulations reported in a later section. Although additional tests on Space Shuttle 
wing leading edge panels, like those reported by Lu at ah, are needed, the cited references 
represent the best data available to the authors at the time the simulations were conducted. 

MATERIAL MODEL 

Composite materials are used in structural , 24 orbital debris shielding , 25 and thermal 
protection 4 applications on a variety of spacecraft, hence their response to hypervelocity 
impact effects has been analyzed in a number of previous experimental 26 and computa- 
tional 27 studies. Previous material modeling work has considered both micromechanical 28 
and anisotropic continuum models. The present paper employs an anisotropic continuum 
approach, since the high computational cost of micromechanical models normally precludes 
their use in structural scale simulations. Large deformation, anisotropic continuum models of 
composite materials 29,30 normally address shock physics problems by extending small strain 
formulations originally developed for applications in structural mechanics. In an alternative 
approach, the present work starts with the finite strain, hybrid particle-element kinematics 
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discussed in an earlier section, and then formulates: (1) an anisotropic strain energy den- 
sity function which depends on a general deviatoric Lagrangian strain tensor, and (2) an 
anisotropic, temperature and rate dependent plastic flow rule which depends on an effective 
stress 31 and satisfies a general isochoric deformation constraint. Both the strain energy den- 
sity function and the plastic flow rule: (1) account for differences in material response under 
tension and compression, (2) account for material reference frame dependence under large 
deformations, and (3) satisfy first and second law thermodynamic constraints. The approach 
used here has been applied with success to model isotropic materials. 32 It is motivated by a 
focus on hypervelocity problems, where large deformation dynamics are of central interest, 
and by the kinematic form of the hybrid numerical method used in the present paper. As 
indicated in the later section on simulations, this material modeling approach provides an 
accurate description of hypervelocity impact effects in reinforced carbon-carbon. Potential 
application of the formulation to model impact in other composites, such as graphite-epoxy 
or Kevlar-epoxy, is of interest for future work. 

In the case of an orthotropic material, with distinct elastic moduli in tension and com- 
pression, the shear strain energy density per unit reference volume is 

= (1 - d) ftd £ fti [(1 + 7i) + (1 - 7i) sgn(Efr )} (ET? + 

2 = 1 
3 3 

(1 - d) ft, £ £ (1 - <%) ft, (EiP) 2 (18) 

4=1 j = 1 

where Sij is the Kronecker delta, /i 0 is a reference elastic modulus, and the parameters 
Hij = fiji are dimensionless constants. The parameters y* are the ratios of the elastic moduli 
in compression to those in tension, while the Ef™ are the components of the elastic shearing 
strain, expressed in a material reference frame. Note that this function is analytic, since 
a change in modulus from tension to compression occurs when the corresponding material 
strain component is zero. 

A plastic flow rule for an anisotropic, rate dependent material, which satisfies the afore- 
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mentioned isochoric plastic deformation constraint, may be obtained by extending a large 
strain Lagrangian formulation previously developed for use in hypervelocity impact applica- 
tions. 32,33 The flow rule is 

E,, = p^fSiM PSP ( 19 ) 

where A is a positive proportionality coefficient, S p is the effective stress, 


S" = m pT m t n t n pT s 


( 20 ) 


and S is the deviatoric stress tensor 

q _ 

<9E e 

The first two coefficients in the ffow rule impose the isochoric plastic deformation constraint, 
and are defined by 

= T= nT^ (c '" r + TC " ) (22) 

and 

gT = T— -tr(T)I (23) 

3 

for any symmetric second order tensor T. The third coefficient performs a component 
transformation from a material reference frame to a fixed global frame, and is defined by 

M T T = R t T R (24) 

for any symmetric second order tensor T. The last coefficient in the ffow rule defines an 
effective stress transformation, in a material reference frame, using 

M pT p = Q (25) 



for symmetric second order tensors P and Q, with component forms 

2 ota Pa 


Qii 


(1 + Pa) + (1 - Pu) sgn(Ef™ ) 


i = 1,2,3 


(26) 
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and 


Qij 


2 OLa Pj 


IJ * IJ 


(1 + fa) + (1 - fa) sgn(J - 1) 


3 


(27) 


The parameter a,j = atji is the ratio of a reference yield stress to the yield stress for the ijth 
stress component, while the parameter flij = (3ji is the ratio of the strength in compression 
to that in tension for the ijth stress component. 

The rate dependent, strain hardening, thermal softening yield function is 


/ 



(28) 


where Y is the yield stress 


Y = l (l-d)Y 0 (1 -k 9 h ) (1 + V £ p ) r 


i + c % ( — 

£os J 


(29) 


with Y 0 the reference yield stress, 77 a strain hardening coefficient, n a strain hardening 
exponent, ( a strain rate hardening coefficient, m a strain rate hardening exponent, e Q a 
reference strain rate, n a thermal softening coefficient, and 6 H the homologous temperature 

e-e n 


e H = ^^Y 

9 m ~ 9 0 


(30) 


where 9 0 and 9 m are reference and melt temperatures. 

In a numerical implementation, the aforementioned plastic flow rule is expressed in in- 
cremental form. That is the incremental plastic strain at each time step is computed using 
the incremental proportionality coefficient 

( ||S P || — Y \ 

AA = max 0 , 

V (1 — d) 2 Ho ) 

The shear damage variable (d) models the transition from an intact to a failed medium, 
evolving from an initial value of 0 to a final value of 1 over a fixed number of time steps 34 when 
any stipulated element failure criterion is satisfied. The simulations discussed in the next 
section incorporate accumulated plastic strain, melt temperature, and maximum compression 
failure criteria, although other criteria may be accommodated. 
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IMPACT SIMULATIONS 


The material model just described was applied in a series of three dimensional simulations 
of hypervelocity impacts on reinforced carbon-carbon. The simulations employed a hybrid 
particle-finite element method and the material properties listed in Tables 1 and 2. An 
initial set of simulations was used to validate the material model, compare results obtained 
using analytic (Mie Gruneisen) and tabular 35 (SESAME 3715) equations of state, and check 
numerical convergence of the simulation results. A second series of simulations was then 
performed to estimate orbital debris impact effects at velocities beyond the range of current 
experimental methods. 

The first set of eight simulations modeled NASA JSC experiments B1028 and B1040, 9 
which involved oblique impacts of aluminum spheres on reinforced carbon-carbon target 
plates at a velocity of seven kilometers per second. The target plates were 0.63 cm in 
thickness, including upper and lower surface coatings composed of silicon carbide, each 
0.08 cm in thickness. Table 3 lists the simulation parameters, including projectile diameter 
(D), impact velocity (u), impact obliquity (0, with zero degrees a normal impact), number 
of elements spanning the target thickness (N e ), and the equation of state used to model 
the aluminum projectile. Tabular equation of state data was not available for the target 
materials. 

Figures 1 through 4 show example plots for a simulation of experiment B1028. Figure 
1 shows the initial configuration while Figures 2 through 4 show the simulation results at 
50 microseconds after impact. The sectioned plot in Figure 4 depicts plate perforation 
and coating spall similar to that observed in the corresponding experiment. Table 3 lists 
simulation results for the diameter of the RCC perforation (D p ) and the average diameter 
of the target region over which the silicon carbide coating was removed (D c ). The results of 
the validation simulations suggest the following conclusions: 
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• the material model developed here can provide good estimates of both the RCC per- 
foration diameter and the extent of the spalled coating region, for oblique impacts at 
seven kilometers per second, 

• accurate estimates of the RCC perforation diameter require a mesh resolution sufficient 
to place 8 elements across the target plate, 

• accurate estimates of the diameter of the spalled coating region require a mesh resolu- 
tion sufficient to place 16 elements across the target plate, and 

• the simulation results are not sensitive to the choice of projectile equation of state. 

Table 4 shows the relative computational cost of simulations of experiments B1028 and 
B1040 run at three different mesh densities. As is well known, in three dimensional models 
the particle count increases with the cube of the increase in resolution, while the time step 
decreases linearly with the increase in resolution, so that the total computational cost of 
high resolution models is considerable. 

A second set of twelve simulations was performed to investigate orbital debris impact 
effects at velocities beyond the current experimental range. The simulations involved spher- 
ical aluminum projectiles, at three different projectile diameters, an impact obliquity of 30 
degrees and impact velocities of 7, 10, and 13 km/s. In the case of the largest projectile, 
simulations were performed using both an analytic and a tabular equation of state. The 
target assumed in these simulations was identical to that involved in the aforementioned ex- 
periments. In the target mesh 8 elements spanned the plate thickness, so that the resolution 
level was sufficient to estimate the diameter of the RCC perforations, but not the extent of 
the region of coating spall. Figure 5 shows simulation results for the diameters of the RCC 
perforations, as a function of projectile diameter, impact velocity, and projectile equation 
of state. In Figure 5, MG denotes the Mie Gruneisen analytic equation of state, while SES 
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denotes the SESAME tabular equation of state. The results of these simulations suggest the 
following conclusions, for the impact velocity range and impact obliquity considered: 

• perforation diameters increase with both projectile size and impact velocity, over the 
full range of the simulations, 

• for a fixed projectile size, perforation diameters increase with impact velocity at an 
approximately linear rate, 

• for a fixed impact velocity, perforation diameters increase with projectile size, but at 
a declining rate, and 

• the simulation results are not sensitive to the choice of projectile equation of state. 

Note that Figure 5 is not a ballistic limit plot; rather it plots perforation diameter versus 
impact velocity, so that the indicated trends are not unexpected. 

Although the preceding results are informative, they consider only a limited rage of 
projectile size and obliquity. Hence the scaling of the simulation results, as compared to the 
available experimental data, is of considerable interest. Figure 6 shows a plot of perforation 
diameter versus normal impact momentum for the 11 different projectile size and impact 
velocity combinations modeled in the present computational study, as well as corresponding 
data for 15 published experiments. 9,36 The experiments involved projectile diameters ranging 
from 0.039 to 0.628 cm, impact velocities ranging from 2.49 to 7.32 km/s, and impact 
obliquities ranging from 0 to 80 degrees. The simulations involved a more limited range 
of projectile diameters (0.123 to 0.360 cm) and obliquities (30 to 45 degrees), but a much 
higher range of impact velocities (7 to 13 km/s). All of the simulations and experiments 
of course involved the same target configuration. The data in Figure 6 suggests that the 
experimental and simulation results for the diameters of RCC perforations scale with normal 
impact momentum in a similar fashion. Although these results do not establish a universal 
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scaling relation for the problem of interest, they do suggest that the scaling of perforation size 
with normal impact momentum observed in experiments below 8 km/s may be extrapolated 
to much of the velocity range of interest in orbital debris impact applications. 

CONCLUSION 


The present paper has formulated an anisotropic, rate dependent material model for use in 
the simulation of hypervelocity impact problems. The material model was developed to study 
orbital debris impact effects on reinforced carbon-carbon materials, and has been validated in 
simulations of hypervelocity impact experiments conducted at 7 km/s. The validated model 
was applied to simulate impacts at velocities beyond the experimental range. The results 
indicate that momentum scaling analysis, used to correlate a wide range of experiments 
below 8 km/s, has application in predicting perforation diameters for reinforced carbon- 
carbon targets at velocities as high as 13 km/s. The ability of reinforced carbon-carbon 
to retain its strength at high temperatures suggests that accurate strength models of this 
material are important in simulations of impact effects over the entire orbital debris velocity 
range. 

Some conclusions relevant to future work are suggested: 

• additional high resolution simulations are needed in order to investigate the spallation 
of silicon carbide coating at velocities above the current experimental range, 

• additional mechanical properties testing is needed, at elevated temperatures and high 
strain rates, to support the development and validation of improved strength models 
for reinforced carbon-carbon, 

• additional equation of state research is needed, to provide tabular data applicable to 
reinforced carbon-carbon materials over a wide range of impact velocities, and 
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• the development of advanced thermal protection materials should in the future include 
experimental work aimed at detailed characterization of their mechanical as well as 
their thermal properties. 

The use of composite materials in spacecraft applications complicates both experimental 
and computational studies of impact effects. High cost materials with long fabrication lead 
times, such as reinforced carbon-carbon, and the limitations of current experimental impact 
techniques motivate the increased use of computer simulation in the design of spacecraft for 
micrometeoroid and orbital debris impact effects. 
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Table 1. Material properties 


Material property 

Aluminum 

Silicon Carbide 

RCC 

Reference bulk modulus ( Mbar ) 

0.784 

2.21 

0.0576 

Reference shear modulus (Mbar) 

0.271 

0.240 

0.0718 

Reference soundspeed ( cm /nsec -1 ) 

0.539 

0.829 

0.191 

Mie-Gruneisen gamma 

1.97 

0.95 

0.24 

Mie-Gruneisen slope 

1.34 

1.21 

1.33 

Reference density ( g cm~ 3 ) 

2.70 

3.21 

1.58 

Reference yield stress ( kbar ) 

2.90 

0.771 

0.771 

Specific heat ( bar cm 3 g~ x °K~ 1 ) 

8.84 

7.12 

7.12 

Strain hardening coefficient 

125 

10 

2 

Strain hardening exponent 

0.1 

1.0 

1.0 

Strain rate hardening coefficient 

0.0 

0.0 

0.1 

Strain rate hardening exponent 

0.0 

0.0 

1.0 

Reference strain rate (sec -1 ) 

0 

0 

0.01 

Thermal softening coefficient 

0.567 

0.0 

-1.0 

Melt temperature (°K) 

1,220 

3,840 

3,840 

Maximum compression 

100 

100 

100 

Plastic failure strain 

1.00 

0.10 

0.50 


Table 2. Material model parameters 


Parameters 

Silicon Carbide 

Reinforced carbon-carbon 

II 

bo 

II 

CO 

0.10 

1.00 

II 

to 

to 

II 

-£ 

CO 

CO 

II 

CO 

1.00 

1.00 

/R2 = /^23 

1.00 

0.50 

an = 022 = «33 = «13 

1.00 

1.00 

ai2 = Oi 23 

1.00 

3.73 

/Til = 7?22 = /?33 

2.00 

2.00 

/3l2 = /3l3 = /3 23 

2.00 

2.00 
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Table 3. Numerical results, simulations of NASA JSC experiments B1028 and B1040 


Test 

number 

D 

(cm) 

V 

(km/s) 

0 

(deg) 

N e 

Equation 
of state 

D p 

(cm) 

Error 

(%) 

D c 

(cm) 

Error 

(%) 

B1028 

0.628 

7.01 

45 

8 

Mie Gruneisen 

2.60 

10.3 

3.74 

15.0 





8 

SESAME 3715 

2.65 

8.6 

3.60 

18.2 





16 

Mie Gruneisen 

2.66 

8.3 

4.05 

8.0 





24 

Mie Gruneisen 

2.67 

7.9 

4.08 

7.3 

B1040 

0.478 

6.96 

30 

8 

Mie Gruneisen 

2.12 

3.6 

2.95 

21.3 





8 

SESAME 3715 

1.97 

10.5 

2.95 

21.3 





16 

Mie Gruneisen 

2.00 

10.0 

3.38 

9.9 





24 

Mie Gruneisen 

2.10 

4.5 

3.48 

7.2 


Table 4. Computer resource requirements, simulations of NASA JSC experiments 


Test 

number 

N e 

Total 

particles 

(million) 

Total 

elements 

(million) 

Number of 
processors 

Wall clock 
hours 

B1028 

8 

0.078 

0.036 

16 

14 


16 

0.572 

0.275 

64 

65 


24 

1.857 

0.905 

64 

340 

B1040 

8 

0.078 

0.036 

32 

5 


16 

0.572 

0.275 

64 

74 


24 

1.856 

0.905 

64 

347 
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Figure 1. Initial configuration, simulation of NASA JSC experiment B1028. 
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Figure 2. Particle-element plot of the simulation results at 50 microseconds after impact. 


160 


Of e 



Figure 3. Element plot of the simulation results at 50 microseconds after impact. 
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Figure 4. Sectioned element plot of the simulation results at 50 microseconds after impact. 
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Figure 5. Simulation results for perforation diameter versus impact velocity. 
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Figure 6. Perforation diameter versus normal impact momentum, for hypervelocity impact 

in reinforced carbon-carbon. 
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